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Q . Using the QCD Operator Product Expansion, we derive the real part of the transverse 

^ ! and longitudinal vector vector correlation function with the p, oj quantum numbers to 

QQ \ leading order in density and in at —u"^ —>■ oo. To dimension 6, only twist-2 and 4 

^ ' operators contribute. These OPE, through the energy dispersion relation, provide model 

^ ■ independent constraints for the q dependence of the vector meson spectral density in 

I . nuclear medium. We further make a QCD sum rule type of analysis to extract the mo- 

y I mentum dependence of the vector meson dispersion relation in medium. The contributions 

pH I from twist-2 operators are added up to infinite order to check the validity of the OPE at 

' the relevant Borel window. 

X ■ 

' 1- Introduction 



The properties of vector meson in nuclear medium have been the focus of current in- 
terest due to its potential role to provide one with a direct observable of the nuclear 
medium effects, associated with chiral symmetry restoration, via dileptons in H-A or A-A 
reactions!^. Indeed dileptons from Relativistic Heavy Ion Collisions (RHIC)0 seemed 
to suggest a non-trivial change of the vector meson spectral density in a hot/dense envi- 
ronment, which can be understood in terms of model calculations based on decreasing 
vector meson masses in hot/dense medium|Q,|^. However, model calculations based 
on changes of the vector meson spectral densities using hadronic variables also seem to ex- 
plain the main features of the CERES data. In all of the approaches, the central question 
is, how the spectral density changes in hot/dense matter p]. In this talk, I will provide 
constraints on the three momentum (q) dependence of the vector meson spectral density 
in nuclear medium [p|,[lO|. The approach is based on the Operator Product Expansion 
(OPE) in QCD and provides model independent constraints that can be used to check 
the validity of any model calculation. In the second stage, I will make a simple ansatz for 
the spectral density near the vector meson mass and look at what the parameters of the 
ansatz should be to satisfy the constraints 0. 



AvH Fellow. 
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2. OPE 

Consider the correlation function of the vector current = q'jfj^q in nuclear matter; 

n^.(u;,q) = zJd'xe''^^G\T[J,{x)MmG)- (1) 

Here \G) is the nuclear ground state at rest and q = {uj,q). In what follows, when we 
give result for explicit vector meson, we will use the currents J^'*^ = |('U7^m =f d'y^d) for 
the p, uj mesons. 

In general, because the vector current is conserved, the polarization tenser in eq. (|ID will 
have only two invariant functions [ITT |. 



Ii,,{uj, q) = lirq^Vl, + n^g'?;,, (2) 



00 ~ Oi ~ iO 



where we assume the ground state to be at rest, such that, Pqq = Poi — 
^Ij = ^ij - q^qj/q^ and P^^ = {q^,qy/q^ - g^,u - P^^). When q ^ 0, 11^ = IIt, as in the 
vacuum. 

In this work, we are only interested in the small three momentum (q) dependence of 
the polarization functions. Therefore we will make a small q expansion of the correlation 
function and look at its energy dispersion relation at fixed q, 

Re^L,T(c^^ q^) = Re (n°(cj^ 0) + ^^,^(^^ 0) q^ + 

T~2 — + — q + ■■ ' 3) 

yiw^ — uj^j [u^ — uj^j J 

where p{u,q) = l/7rImn-^(M^, q^), and R denotes the retarded correlation function. The 
OPE for n° and its change in nuclear medium together with its corresponding energy 
dispersion relation provide constraints on the mass shift and spectral changes at q = 0[^]. 
In this talk, we will discuss the three momentum dependence by studying the OPE and 
its changes in medium for Hit- 

In general, the OPE |]T2|,[TB| for the polarization function at = —lJ^ + q^ — > large 
will look as follows. 



n^;.(t^,q) = {q^iqu- 9i,vq^ 



n ^ 



+ ^[-9^luq^llq^^2 + S-mmi 9^9/^2 + q^^q^il9u^,2 + 9^i^ii9vti2Q'^ 

T = 2 fc = l 

r)2k 

n ■ -n A2k+r,T 



/q']q,s 



■)2fc 



T=2 k=l 



Q 



Q2k+7 



(4) 



Here we have pulled out the trivial ^ dependence so that A'^''^, C"^'"^ represents the 
residual Wilson coefficient times matrix element of an operator of dimension d and twist 
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T = d — s, where s = 2k is the number of spin index of the operator. The first set of terms 
comes from the OPE of scalar operators, the second set from operators with spin, which 
have been written as a double sum of twist and spin. A, C represents the two linearly 
independent sum of operators, which reflect the two linearly independent polarization 
directions. 

The longitudinal and transverse polarization functions can be obtained from eq.(|^). 
The q dependence coming from the first line of eq.(^), namely the contribution from the 
scalar operators, comes from the q dependence in Q^. This forms the so called "trivial" 
q dependence, and comes from replacing ^ u'^ — q^ when going from zero to finite 
three momentum. Here, we are not interested in these trivial dependence and will only 
look at the "non-trivial" q dependence in Il\ j,. Consequently the scalar operators do not 
contribute to Il\ j,. Only operators with spins contribute to the non-trivial q dependence. 
A prescription to find the nontrivial q^ dependence in the OPE is to first calculate the total 
contribution proportional to q^ and then subtract out the trivial dependence; ^{d)^ 

{d — f)^;;^, and its contribution to 11^ would be {d — |)-;r+2. Another equivalent method 
is to express the polarization function in terms of n(Q^, q^) and extract the linear q^ 
term. 



3. Linear density Approximation 

In the linear density approximation, 

{G\A\G) = m\^) + ^{p\m. (5) 

where the first term denotes the vacuum expectation value, which vanishes for operators 
with spin, and the second term the nucleon expectation value with the normalization 
{p\p) = 2po'5^(p — p)- Pn,^ denotes the nuclear density and the mass of the nucleon 
respectively. 

As in the vacuum, we will truncate our OPE at dimension 6 operators. This implies 
that in our OPE in eq.(D, we will have contributions from (r, s) = (2, 2), (2, 4), (4, 2). 
The nucleon matrix elements of the r = 2 operators are very precisely known. The r = 4 
matrix elements appearing in the p, u sum rule are similar to those appearing in electron 
DIS[n] and have been estimated |IH|TB|up to about ± 30% uncertainty from available 



DIS data from CERN and Slac. 

The final form of the OPE for 11^ looks as follows. 

ni,.M/Pn = ^ + ;^. (6) 
For p,u, the transverse (T) and longitudinal (L) parts give, 

+\rn (-(1 + (3){K' + ^-K' + ^K^) + K^l ± 1) 
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62 = -Ici^mAr' - C^.mAl (7) 

63 = ilcl,-lcl,)m'Al-^' + {C^^,-5Cl,)m'A^ 

+ (8) 
where ± refers to the p and u case. Here, for even n, 

Al = 2 [\xx"-'[q{x,Q^) + qix,Q^)] , A^^ = 2 dxx''~^G{x,Q^), (9) 

JO JO 

where Q^) and (^(a;, Q^) are the quark and gluon distribution functions. We will 
use the HO parameterization for these obtained in reffl^ which should be used with 
the Wilson coefficients C's in the MS scheme [p!8[] and given below. Terms proportional 



Table 1 

Wilson coefficients for r = 2. T{R) = |, C2(i?) = |. 

Wilson coefficient n = 2 n = 4 

~Cl^^TT^B^ 1 + ^(0.44) l + t(6.07) 

(Jl _ (fl) 4 as 16 as 16 



r^G — as_T( T^\^\—^ ^ _|_ J nf+n+2_/i _|_ v^n iNl as_( _]_\ as ( 133 \ 

'-"2,n iiT-^ f\-{n+l) (n+2) n2 n(n+l){n+2) V ^j=l j /J 47r V 2^ A-k^ im' 

r^G _ as '-pl TD\ 1 16 2 as 4 

'-^L,n 47r / (n+l)(ra+2) 47i_3 Air 15 



to K' s come from r = 4, s = 2. We will use the set of K values obtained in ref.[p!5 
( K\K^,Kl^,K9) = (-0.173GeV^0.203GeV^-0.083GeV^-0.238GeV^) and we take 
(3 = 0.5. For the r = 4 operators, we neglect the dependence. 

4. Constraints 

We will now make a Borel transformation of u"^ times eq. (^) and look at the dispersion 
relation. 

j dsp\s)se''"'" +hscaU. (10) 




M2 2M4 



The reason for multiplying by uP is to get rid of any possible subtraction constants propor- 
tional to 1/a;^. There could be still further subtraction constants proportional to hscatt/'^^-i 
which we included, hscatt should also be calculated in any model calculations. The value 
from particle- hole intermediate state gives hscatt = — l/(4^)Pn for the longitudinal p,uj 
meson and zero for the transverse parts. 

After the Borel transformation, the contribution from dimension=?T, operators in the 
OPE side is further divided by (n — 2)!. This will make our truncation at dimension 6 
operators valid even to smaller Borel mass region. In fig.l, we have plotted OPE as a func- 
tion of M^. The solid line shows the left hand side of eq. (|lO|) . The dot-dashed line without 
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the twist-4 contribution and the long dashed hne without dimension 6 contribution. At 
^min ~ -SGeV^ (^mm ~ 2.5GeV^), the contribution of dimension 6 operators become 
about 40% of the contributions from dimension 4 operators for transverse (longitudinal) 
polarization. Therefore, the OPE should be reliable only above these minimum values 
of the Borel mass, where the contributions from higher dimensional operators should be 
similarly suppressed. To support our argument, we have also plotted (short dashed line) 
the contributions from adding up the infinite sum of the twist-2 operators. This is pos- 
sible because for the twist-2 operators, all the moments are know. As can be seen from 
the figure, once we are above the minimum Borel mass, the difference to the solid line is 
vanishingly small, confirming the validity of our OPE above the Minimum Borel mass. 
Now the constraints for the spectral density would be eq.(p!0|) , but applied only above 




Figure 1. Solid line denotes our OPE, short dashed line with the infinite sum of the 
twist-2 operators, dashed line OPE without dim 6, dot-dashed OPE without twist-4. All 
curves are at nuclear matter density and GeV^ unit for both axis 



mm 

As can be seen in eg . ([I0|) , the Borel transformation also changes the weighting factor 
of the spectral density to an exp{—s/M'^). This has the following advantage for practical 
applications of our constraint. For small values of the Borel mass, the contribution of the 
spectral density at larger energy is exponentially suppressed. Consequently, in a model 
calculation, one can concentrate on the change of the spectral density near the vector 
meson mass region and below and model the higher energy changes with a simple pole 
like contribution. 

4.1. Results with simple ansatz 

We now show the result with the following ansatz of the spectral density. 

FTpinip + ^q2) 

p{u) = — i'"%^r2 I ^^+c6{u^-{So + sq^)). (11) 
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We expand this to linear order in to obtain and use the constraints in eq.(0) to 
obtain the best fit value of the constants f,a,g,s. In the limit F — > and g = 0, the 
results were obtained in ref. 0]0 which translates into the following momentum dependence 
in the mass. ^^^^ = 1 - (0.023 ± 0.007)(^)2(^), where q is in GeV unit and po is the 
nuclear matter density. For the u meson, 0.023 changes to 0.016. 

5. Conclusion 

We have shown that there are model independent constraints on the three momentum 
dependence of the vector meson spectral density in nuclear medium. Any model calcula- 
tions can be checked to see if it satisfies the constraints. Phenomenological models based 



on p-wave nucleon resonances seems to overestimate the momentum dependence |T9 
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